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Abstract. A model for Gamma ray bursts inner engine based on quark stars (speculated to exist in nature) is 
presented. We describe how and why these objects might constitute new candidates for GRB inner engines. At 
the heart of the model is the onset of exotic phases of quark matter at the surface of such stars, in particular the 
2-flavor color superconductivity. A novel feature of such a phase is the generation of particles which are unstable 
to photon decay providing a natural mechanism for a fireball generation; an approach which is fundamentally 
different from models where the fireball is generated during collapse or conversion of neutron star to quark star 
processes. The model is capable of reproducing crucial features of Gamma ray bursts, such as the episodic activity 
of the engine (multiple and random shell emission) and the two distinct categories of the bursts (two regimes are 
isolated in the model with ~ 2 s and ~ 81 s burst total duration). 
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1. Introduction 

A central problem contributing to the Gamma-ray bursts 
(GRBs) mystery is the unknown nature of the engine pow- 
ering them (Kouveliotou et al. 1995; Kulkarni et al. 1999; 
Piran 1999a; Piran 1999b). Many have been suggested but 
it is fair to say that we are still far from a definite an- 
swer. Regardless of the nature of the engine, however, it is 
widely accepted that the most conventional interpretation 
of the observed GRBs result from the conversion of the ki- 
netic energy of ultra-relativistic particles to radiation in 
an optically thin region. The particles being accelerated 
by a fireball mechanism (or explosion of radiation) taking 
place near the central engine (Goodman 1986; Shemi & 
Piran 1990; Paczyhski 1990). 

The first challenge is to conceive of circumstances that 
would create a sufficiently energetic fireball. Conversion of 
neutron stars to quark stars (Olinto 1987; Cheng & Dai 
1996; Bombaci & Datta 2000) has been suggested as one 
possibility. Other models also involve the compact object 
element; such as black holes (Blandford & Znajek 1977) 
and coalescing neutron stars (Eichler et al. 1989; Ruffert 
& Janka 1999; Janka et al. 1999). We show in this work 
that the plausible existence of quark stars combined with 
the onset of a newly revived state of quark matter - called 
color superconductivity - in these objects offers a new way 
of tackling the GRB puzzle (Ouyed 2002). Here we will 
argue that quark stars might constitute new candidates 
for GRB inner engines. 
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Quark matter at very high density is expected to be- 
have as a color superconductor (see Rajagopal & Wilczek 
2000 for a review). Associated with superconductivity is 
the so-called gap energy A inducing the quark-quark pair- 
ing and the critical temperature (Tc) above which thermal 
fluctuations will wash out the superconductive state. A 
novel feature of such a phase is the generation of glueball 
like particles (hadrons made of gluons) which as demon- 
strated in Ouyed & Sannino (2001) immediately decay 
into photons. If color superconductivity sets in at the sur- 
face of a quark star the glueball decay becomes a natural 
mechanism for a fireball generation. 

The paper is presented as follows: In Sect. 2 we briefly 
describe the concept of color superconductivity in quark 
matter. Glueball formation and their subsequent two- 
photon decay is described. Sect. 3 deals with quark stars 
and the onset of color superconductivity at their surface. 
In Sect. 4, we explain how GRBs are powered in this pic- 
ture and show that variability (multiple shell emission) is 
inherent to the inner engine. We isolate two GRB regimes 
in Sect. 5 associated with small and massive quark stars. 
The model's features and its predictions are summarized 
in Sect. 6 while a discussion and conclusion follows in Sect. 
7 where the model's assumptions and limitations are high- 
lighted. 



2. Color Superconductivity 

While in this paper we deal mostly with the astrophysics 
aspect of the model, we nevertheless give a brief overview 
of color superconductivity and the glueball-to-photon de- 
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Fig. 1. A schematic representation of a possible QCD 
phase diagram (Rajagopal & Wilczek 2000). At high tem- 
perature and density, matter is beheved to be in a quark- 
gluon plasma phase (QGP). The hadronic phase lies in 
the region of low temperature and density. At very high 
density but low temperature, when nuclei melt into each 
other, it has been suggested that a color superconductive 
phase might set in. 2SC denotes a 2-flavor color super- 
conductive regime. The arrow depicts a plausible cooling 
path of a HQS surface leading to the onset of color super- 
conductivity. 



cay process which leads to the fireball. The interested 
reader is referred to Ouyed & Sannino (2001) for the un- 
derlying physics. For a recent review see Sannino (2002) 

2.1. 2-flavour color superconductivity 

A reasonable Quantum Chromo-Dynamics (QCD) phase 
diagram (in the fJ-—T plane, where fi is the chemical poten- 
tial simply related to matter density) is shown in Fig. |l|. 
At high temperature and density, matter is believed to 
be in a quark-gluon plasma phase (QGP). The hadronic 
phase lies in the region of low temperature and density. At 
high densities but low temperatures, when nuclei melt into 
each other, it is now believed that a color superconductive 
phase sets in. This phase is characterized by the formation 
of quark-quark condensate. In the 2-flavor color supercon- 
ductivity (2SC) the up and down quark come into play 
during pairing. Furthermore, 2SC is characterized by five 
out of the eight gluons acquiring mass. We refer the inter- 
ested reader to Rajagopal & Wilczek (2000) for a review 
of the dynamical properties of 2SC. 

2.2. Light GlueBalls 

The 3 massless gluons in the 2SC phase which bind into 
light glueballs (LGBs) together with the quarks up and 
down constitute the 2SC phase mixture. In Ouyed & 
Sannino (2001) we studied certain properties of these 
LGBs. Among the properties relevant to our present study 
we found, i) The LGBS decay into photons with an asso- 
ciated lifetime of the order of 10"^** s; ii) The mass of the 
LGBs is of the order of 1 MeV. 



3. Quark stars 

We now turn to study the astrophysical consequences 
when such a state sets in at the surface of quark stars. 
As such, we first assume th at q uark stars exists in nature 
(further discussed in Sect. 7T) and constitutes the first 
major assumption in our model. 

3.1. Hot Quark stars 

We are concerned with quark stars born with surface tem- 
peratures above Tc. We shall refer to these stars as "hot" 
quark stars (HQSs) in order to avoid any confusion with 
strange stars which are conjectured to exist even at zero 
pressure if strange matter is the absolute ground state of 
strong interacting matter rather than iron (Bodmer 1971; 
Witten 1984; Haensel et al. 1986; Alcock et al. 1986; Dey 
et al. 1998). 

We borrow the language of the MIT-bag model for- 
malism at low temperature and high density to describe 
HQSs (Farhi & Jaffe 1984). This gives a simple equation 
of state 



P = 6(p-pHQs)c^, 



(1) 



where 6 is a constant of model-dependent value (close to, 
but generally not equal, to 1/3 of the MIT-bag model), 
and phqs is the density at zero-pressure (the star's surface 
density). We should keep in mind that Tc/fi <C 1 as is 
confirmed later. 

Features of HQSs are - to a leading order in Tc - iden- 
tical to that of strange stars. The latter have been studied 
in details (Alcock et al. 1986; Glendenning & Weber 1992; 
Glendenning 1997). Of importance to our model: 

i) The "surface" of a HQS is very different from the sur- 
face of a neutron star, or any other type of stars. Because 
it is bound by the strong force, the density at the surface 
changes abruptly from zero to phqs- The abrupt change 
(the thickness of the quark surface) occurs within about 
1 fm, which is a typical strong interaction length scale. 

ii) The electrons being bound to the quark matter by 
the electro-magnetic interaction and not by the strong 
force, are able to move freely across the quark surface 
extending up to ~ 10'^ fm above the surface of the star. 
Associated with this electron layer is a strong electric field 
(5 X 10^"^ V/cm)- higher than the critical value (1.3 x 10^^ 
V/cm) to make the vacuum region unstable to sponta- 
neously create (e+,e~) pairs. 

iii) The presence of normal matter (a crust made of 
ions) at the surface of the quark star is subject to the enor- 
mous electric dipole. The strong positive Coulomb barrier 
prevents atomic nuclei bound in the nuclear crust from 
coming into direct contact with the quark core. The crust 
is suspended above the vacuum region. 

iv) One can show that the maximum mass of the crust 
cannot exceed Mcrust — 5 x 1O~^M0 set by the require- 
ment that if the density in the inner crust is above the neu- 
tron drip density (pdrip — 4.3 x 10^^ g/cc), free neutrons 
will gravitate to the surface of the HQS and be converted 
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to quark matter. This is due to the fact that neutrons 
can easily penetrate the Coulomb barrier and are readily 
absorbed. 



3.2. Cooling and 2SC layer formation 

The HQS surface layer might enter the 2SC phase as illus- 
trated in Fig. |l|. In the QCD phase diagram (Fig. ||), (pBo i 
Tbo) is the critical point beyond which one re-enters the 
QGP phase (the extent of the 2SC layer into the star) . The 
star consists of a QGP phase surrounded by a 2SC layer 
where the photons (from the LGB/photon decay) leaking 
from the surface of the star provides the dominant cooling 
source. This picture, as illustrated in Fig. ^ is only valid if 
neutrino cooling in the 2SG phase is heavily suppressed as 
to become slower than the photon cooling. Unfortunately, 
the details of neutrino cooling in the 2SC phase is still a 
topic of debate and studies (Carter & Reddy 2000; Schaab 
et al. 2000 to cite only few). One can only assume such 
a scenario which constitutes the second major assump- 



tion in our model. In Sect. 7.2, we discuss the remaining 
alternative when photon cooling is dwarfed by neutrino 
cooling. 

3.3. LGBs decay and photon thermalization 

The photons from LGB decay are generated at energy 
< Tc and find themselves immersed in a degener- 
ate quark gas. They quickly gain energy via the inverse 
Compton process and become thermalized to Tc- We es- 
timate the photon mean free path to be smaller than few 
hundred Fermi (Rybicki & Lightman 1979; Longair 1992) 
while the 2SC layer is measured in meters (see Sect. 5.2). 
A local thermodynamic equilibrium is thus reached with 
the photon luminosity given by that of a black body radi- 
ation. 



= 3.23 X 10 



52 



ergs s 



5 km' '10 MeV 
The energy for a single 2SC event is thus 



(2) 



(3) 



where M2SC = '^2SC-^^hqs is the portion of the star in 
2SC. Here, 62SC depends on the star's mass while (5lgb 
represent the portion of the 2SC that is in LGBs (intrinsic 
property of 2SC; see Ouyed & Sannino 2001). The emis- 
sion/cooling time is then 



Atcool — 



with 



32SCOLGB- 



(4) 



4. Powering Gamma-Ray Bursts 

4.1. Fireball and baryon loading 

The fireball stems from the LGB decay and photon ther- 
malization. The photons are emitted from the star's sur- 
face into the vacuum region beneath the inner crust ('^ 10^ 
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Fig. 2. The episodic emission as illustrated in the QCD 
phase diagram. The 2SC front spreads deep inside the 
star and stops at Bq before re-entering the QGP phase. 
Following photon cooling, heat flows from the core and 
re-heats the surface. The star then starts cooling until 
Al is reached at which point the stage is set for the 
2SC/LGB/photon process to start all over again (Ai 
Bi) resulting in another emission. 

fm in size). Photon-photon interaction occurs in a much 
longer time than the vacuum region crossing time. Also, 
the cross-section for the creation of pairs through interac- 
tions with the electrons in the vacuum region is negligible 
(Rybicki & Lightman 1979; Longair 1992). The firebaU 
energy is thus directly deposited in the crust. If its en- 
ergy density, aT^ (with a being the radiation density con- 
stant), exceeds that of the gravitational energy density in 
the crust, energy outflow in the form of ions occurs. One 
can show that the condition 



4 GA/hqs 

fl-ic ^ Pcrust 

^HQS 



(5) 



where Pciust is the crust density and G the gravitational 
constant, is equivalent to 



Tr 



30 Mev 



^4 ^ ^^HQS -J ^ 5 km ^ ^ Pcrust 



Mr 



' -RhQS Pdrip 



(6) 



which is always true if Tc > 30 MeV. The fireball is thus 
loaded with nuclei present in the crust. More specifically, 
it is the energy transfer from photons to electrons which 
drag the positively charged nuclei in the process. Note that 
the 2SC layer is not carried out during the two-photon de- 
cay process because of the star's high gravitational energy 
density: pHQs/Pdrip >> 1. 

4.2. Episodic behavior 

The star's surface pressure is reduced following photon 
emission^. A heat and mass flux is thus triggered from 
the QGP phase to the 2SC layer re-heating (above Tc) and 
destroying the superconductive phase. The entire star is 



^ The pressure gradient in the 2SC layer is Ap oc (8 — 5)r* 
(Farhi & JafFe 1984) where the massless gluons (3 out of 8) 
have been consumed by the LGB/photon process. 
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now in a QGP pliase (5 gluons-^ 8 gluons at the surface) 
and hence the coohng process can start again. This corre- 
sponds to the transition [p(Bo), T(Bo)] [phqs, 7"(Bo)] 
in the QCD phase diagram (thermal adjustment). The 
stage is now set for the 2SC/LGB/photon process to start 
all over again resulting in another emission. For the subse- 
quent emission, however, we expect the system to evolve 
to point Bi generally located at different densities and 
temperatures than Bq (see Fig. The cycle ends after 
TV emissions when p^Bn) ~ Phqs- 

The time it takes to consume most of the star (the glue 
component) by this process is 



HQSt 



engine 



ls(^)(4^)^(^)\(7) 



Mrr 



^7 ^y-^o ' Rhqs 

which is representative of the engine's activity. The above 
assumes quick adjustment of the star following each event, 
but is not necessarily the case for the most massive stars. 

4.3. Multiple shell emission 

The episodic behavior of the star together with the result- 
ing loaded fireball (we call shell) offers a natural mecha- 
nism for multiple shell emission if < 30 MeV. Indeed 
from Eq. ^ a higher value would imply extraction 
of the entire crust in a single emission and no loading of 
the subsequent fireballs. Clearly, Tc < 30 MeV must be 
considered if multiple ejections are to occuiQ. 

The fraction (/) of the crust extracted in a single event 

is, 



AM, 



(8) 



The shell is accelerated with the rest of the fireball con- 
verting most of the radiation energy into bulk kinetic en- 
ergy. The corresponding Lorentz factor we estimate to be, 



r 



shell 



/Mcrust 



(9) 



where we used Eq. (||) and Eq. (g). e and / depend on the 
star's mass and characterize the two emission regimes in 
our model. 



The ratio rsheii.n/rsheii.n-i is a function of e (which de- 
pends on the details of the cooling process and the spread 
of the 2SC front) and / (mostly related to the density 
in the crust). With the two parameters varying from one 
emission to another, the ratio can be randomly greater or 
less than 1. As such, the shells will have random Lorentz 
factors and random energies. Faster shells will catch up 
with slower ones and will collide, converting some of their 
kinetic energy to internal energy. 

5. The two regimes 

When the inner crust density is the neutron drip value, 
one finds a minimum mass star of ~ O.OISMq. For masses 
above this critical value, the corresponding crusts are thin 
and light. They do not exceed few kilometers in thickness. 
Matter at the density of such crusts is a Coulomb lattice 
of iron and nickel all the way from the inner edge to the 
surface of the star (Baym et al. 1971). For masses below 
O.O15M0, the crust can extend up to thousands of kilo- 
meters with densities much below the neutron drip. This 
allows us to identify two distinct emission regimes for a 
given Tc (< 30 MeV). 

5.1. Light stars (Mrqs < 0.015Mq) 

These are objects whose average density is ~ Phqs 
(-^HQS — x^HQS/^HQs)- The 2SC front extends deeper 
inside the star {S2SC ^ 1)- The star can be represented 
by a system close to Aq in Fig. ||. Each of the few emis- 
sions (defined by e) is thus capable of consuming a big 
portion of the star. Furthermore, the entire crust mate- 
rial can be extracted in a few 2SC/LGB/photon cycles 

(Pcrust/Pdrip ^1). 

Using Eq. (Q), the few emissions lead to 

ttot - fraction x iongino 

- fraction x 0.25 s (^^^QS )(U^)2(30 MeV,, 



^0.01 MQ"i?HQS' 



(11) 



4.4. Shell-shell collision 

The Lorentz factor for the 



shell is 



shell,! 



/n-^^crust.n 

^ ^ ^ fn— 1-^ (1 ^n—l)^ 

£n-l fn (1 — .fn-l) 



shell, n — 1 



(10) 



^ Even if Tc turns out to be greater than 30 MeV, in which 
case the entire crust will be blown away (Eq. (^), one can 
imagine mechanisms where crust material is replenished. By 
accretion, for instance, if the HQS is part of a binary. There 
are also geometrical considerations where asymmetric emis- 
sion/ejection can occur due to the rapid rotation of Quark 
Stars; here only a portion of the crust is extracted at a time. 
This aspect of the model requires better knowledge of the con- 
ditions and environments where HQSs are formed. 



where ttot is representative of the observable time which 
takes into account the presence of the crust. 

5.2. Massive stars (Muqs > O.OISM© j 

The surface density of a massive star being that of a light 
star (phqs given by P = in Eq. @)), defines a standard 
unit in our model. In other words, the mass of the 2SC 
layer in a massive star case is 



AAf2SC,m — M2SCJ , 



(12) 



where "m" and 'T' stand for massive and light, respec- 
tively. It implies 



AP2SC,m _ 1, Pi ,3 1,1 km 3 



s^p;;) -3W 



(13) 
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For a typical star of 5 km in radius, we then estimate a 
2SC layer of about 15 meters thick (much larger than the 
photon mean free path thus justifying the local thermal 
equilibrium hypothesis). Equivalently, 



Ml 1 km 3 
JdZ ~ ^5 km^ 



0.01 



(14) 



where e is the average value. This naturally account for 
many events (or N fireballs). The average number of fire- 
balls with which an entire star is consumed is thus 



TV 



1 



100 . 



(15) 



Since most of the crust is at densities close to the neu- 
tron drip value, Eq. (|^) implies that only a tiny part of 
the crust surface material (where Pcrust << Pdrip) can be 
extracted by each of the fireballs. This allows for a con- 
tinuous loading of the fireballs. 

The total observable time in our simplified approach 
is thus, 



''engine 



ls( 



Mhqsw 5 km 2.30 MeV 



-)( 



Mq Rhqs 
We isolated two regimes: 



)^(- 



T, 



-r 



(16) 



(i) Light stars short emissions. 

(ii) Massive stars long emissions. 

It appears, according to BATSE (Burst and Transient 
Source Experiment detector on the COMPTON-GRO 
satellite) , that the bursts can be classified into two distinct 
categories (Kouveliotou et al. 1993): short (< 2 s) bursts 
and long (> 2 s, typically ~ 50 s) bursts. The black body 
behavior (T|5) inherent to our model puts stringent con- 
straints on the value of T^. which best comply with these 
observations. Using Tc ~ 10 MeV, from Eq. (|ll|) and Eq. 
(Iq) we obtain in the star's rest frame 



tot 



g^^^MHQS)(5km^2^ 



for massive stars (suggestive of long GRBs) , and 



tt 



.J£HQ^..lkm.2 
^0.01 Mq'^'Rhqs' 



(17) 



(18) 



for light stars (suggestive of short GRBs). There is a clear 
correlation (almost one to one) between the observed burst 
time and the time at which the source ejected the specific 
shell (see Figure 3 in Kobayashi et al. 1997, for example). 
Note that Tc ~ 10 MeV implies that only a portion of the 
crust is extracted. This is also consistent with our previous 
assumption (Tc < 30 MeV) and subsequent calculations. 

When Te ~ 10 MeV, Eq. §) gives Perust/pdrip ^ 1/16. 
For an appropriate crust density profile (using the equa- 
tion of state given in Baym et al. 1971), from Eq. (^ we 
find / ~ 0.01. This implies (making use of Eq. dlj 



r.hoii = 2 X 105( 



)(^) 



0.01" / 



(19) 



For massive stars ther 
< r.heii < 2 X 10^ . 



(20) 



For light stars, where both e and / are close to 
unity, FshoU — 10^. The shells are also heavier than in 
the massive stars case. We thus expect stronger shocks 
from the shell-shell collision resulting in harder bursts. 
Combined with our previous results, this is suggestive of 

(i) Light stars ^ short and hard bursts. 

(ii) Massive stars => long and soft bursts. 

Eq. and Eq. (|l|) is simply Eq. (@) rescaled to the 
appropriate object size. We separated two regimes due to 
intrinsic differences in the engine and the crust. From the 
engine point of view, massive stars generate many more 
emissions when compared to light ones, and no substantial 
reduction of the engine time is expected because of the 
omni-presence of the crust. Another important difference 
is related to the physics of the multiple re-adjustments 
following each event which is more pronounced for very 
massive stars. The latter among other factors is related to 
e which can vary from one event to another. 

6. Features and predictions 

6.1. GRB energies 

The maximum available energy is when the heaviest HQS 
(AfHQS,max — '^Mq) is entirely consumed. That is, 



-BLGB.max =i 4 X lO^'^ CrgS . 

The corresponding GRB energy is 

-BcRB^max ^ 1-6 X lO^'' CrgS , 



(21) 



(22) 



where we used a fiducial conversion efficiency of 40% (Sect. 
Since A^HQS.min < 0.015Afo we conclude that, 



-BLGB.min < 3 X 10^^ CrgS , 

implying 

-BGRB.min < 1.2 X 10^^ CrgS . 

6.2. GRB total duration 

From Eq. (|l^) and Eq. (^ we have 

ttot 81 S , 

for typical massive stars, and 
<tot ^2s , 



(23) 
(24) 



(25) 



(26) 



^ Note that a relativistic loaded fireball is not necessarily 
achieved since Fshcii < 100 at times. Gaps in the GRBs spectra 
are thus expected according to our model. 
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for typical light stars. 

Our estimate of the duration time for the massive 
star case should be taken as a lower limit. As we have 
said, a complete model should take into account star re- 
adjustments. Nevertheless, we can still account for a wide 
range in GRB duration by an appropriate choice of differ- 
ent values of the mass and radius. 

6.3. Peak duration vs energy 

The peak duration (tp) is related to the time measured by 
a clock on the shell via (Fenimore & Ramirez-Ruiz 1999), 



in 



tshcll 



1 



AX 



shell 



2r. 



hell 



2r. 



hell 



(27) 



where Cs is the speed of the shock front crossing the shell 
leading to the burst. We find 



(28) 



In the above we used Eq. (||) and Afsheii — AMdust — 
47ri?^^ygj AXcrust /Ocrust • The expression on the right is oc 
Atcooi as can be seen from Eq. (^ and in our model is 
constant for a given star. The only parameter which is 
directly linked to the shell dynamics and energetics is Cg. 
Shock physics gives (Ouyed & Pudritz 1993) 



(29) 



where Ei^t is the shell's internal energy gained during the 
collision observed as the peak's energy (Ep). That is. 



tp oc Ep 



0.5 



(30) 



in reasonable agreement with the power law dependence 
extracted from temporal vs energy structure in GRBs (in- 
dex that is between —0.37 and —0.46; Fenimore et al. 
1995). 

6.4. Shell dynamics 

Take a shell of thickness AXcrust to be extracted from the 
crust. The upper surface of the shell is extracted first while 
its lower surface lags behind by (c — Usheii)i — ct/ (2rgjjj,[[) 
{t is the time to eject the entire shell in the star's rest 
frame). Taking into account mass conservation and the 
fact that AXsheii = ZF^j^gijAXcrust, it is straightforward 
to show 



1 



Pshell OC ^2 



(31) 



hell 



Interestingly enough, these are the required conditions (in- 
cluding the result from Eq. (po|)) in the internal shock 
model which lead to the highest (up to 40%) conver- 
sion efficiency and the most desirable temporal structure 
(Kobayashi et al. 1997; Mochkovitch et al. 1995). 



7. Discussion and Conclusion 

7.1. Existence and formation of quark stars 

In the last few years, thanks to the large amount of fresh 
observational data collected by the new generation of X- 
ray and 7-ray satellites, new observations suggest that the 
compact objects associated with the X-ray pulsars, the X- 
ray bursters, particularly the SAX J1808. 4-3658, are good 
quark stars candidates (see Li et al. 1999). While these 
observations/measurements are hints that such objects 
might exist in nature it remains to explain their forma- 
tion. More importantly to our model, the bimodal mass 
distribution remains to be explained. 



7.1.1. 



lassive stars 



For the massive stars the conversion of neutron stars to 
quark stars is one plausible scenario (Cheng & Dai 1996; 
Ma 1996; Bombaci & Datta 2000). They could also form 
via the direct mechanism following a supernova collapse 
where the core collapsed to a stable quark matter instead 
of neutron matter (Gentile et al. 1993; Dai et al. 1995). 
Both mechanisms would lead to the formation of quark 
stars (strange stars to be more specific) with masses in 
the solar mass range. 

7.1.2. Light stars 

The formation of small quark stars has already been 
discussed (early discussions can be found in Alpar 
1987; Glendenning et al. 1995; see also Chapter 10.5 in 
Glendcnning 1997) although these remain less understood 
than the massive ones. In the case of 4U 1728-34 (where 
a mass of much less than I.OMq was derived; Bombaci 
1999), it seems that accretion-induced collapse of white 
dwarfs is a favored formation mechanism. If the quark 
star formed via the direct conversion mechanism then it 
required too much mass (at least ^ O.8Af0 to be ejected 
during the conversion). 

How and why stars in the O.OIM© range would form 
remains to be explained. Our arguments were solely based 
on theoretical considerations related to the critical density 
in the inner crust (neutron drip) as to differentiate be- 
tween small stars with thick and heavy crust versus stars 
with thinner and lighter crusts. 

7.2. Neutrino cooling and HQSs 

If neutrino cooling is shown to remain efficient in the 2SC 
phase (for comparison of cooling paths between quark 
stars and neutron stars and the plausible effects of 2SC 
on cooling we refer the interested reader to Schaab et al. 
1997; Blaschke et al. 2000; Blaschke et al. 2001), we would 
be left with the scenario where the entire HQS enters the 
2SC phase, in which case the 2SC/LGB/photon process 
(the fireball) occurs only once and inside the entire star. 
Here, one must involve more complicated physics (such as 
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that of the crust) to account for the episodic emissions 
so crucial to any model of GRBs. It is not clear at the 
moment how to achieve this and is left as an avenue for 
future research. 

7.3. 2SC-II stars 

The 2SC/LGB/photon process might proceed until one is 
left with an object made entirely of 2SC. We name such 
objects 2SC-II stars^ which are still bound by strong in- 
teractions (their density is constant ^ Phqs)- 2SC-II stars 
carry an Iron/Nickel crust left over from the GRB phase. 
The crust mass range is < M2sc, crust < 5 x IQ^^Mq 
depending on the efficiency of crust extraction/ejection 
during the GRB phase. 

BATSE observes on average one burst per day. This 
corresponds, with the simplest model - assuming no cos- 
mic evolution of the rate - to about once per million years 
in a galaxy (Piran 1999a). In the Milky way we thus ex- 
pect up to 10^ of 2SC-II stars. Nevertheless, they are tiny 
enough (M < IQ-^Mq, R<1 km) to be difhcuh to de- 
tect. 



0.2 



0.15 



0.05 









GRB991208 


DJ .' 

■S ■' 
■3 ■ 

lu ;' 




'o , ' 
^ : 
■■' 

o ; 




^ ;■ 


// - 

A / 

/ 




/ GRB000418" 







12 3 

Radius (Km) 



Fig. 3. The Mass— Radius plane derived in our model us- 
ing few existing GRBs with known energies and total du- 
ration. The sohd curve shows the Mhqs = ^Phqs^hqs 
equation for phqs — 9pN- 



* The "11" in 2SC is a simple reminder of the final state of 
the star, namely the 2SC with only 5 gluons. 



7.4. The Mass- Radius Plane 

Take observed GRBs with known energies and total du- 
ration. From the burst total energy i^cRB — 0.4Mc^ we 
derive the mass while the burst total duration (ttot) gives 
us the radius (using Eq. (|l|) with ~ 10 MeV). In Fig. | 
we plot the resulting Mass— Radius. Note that while neu- 
tron stars, can only exists above a certain mass (~ O.IMq; 
Baym et al. 1971), there is no lower limit to the mass of 
quark stars. These would be bound by the strong interac- 
tion even in the absence of gravity. 

The solid curve shows the Mhqs = itPhqs-Rhqs equa- 
tion which is a reasonable approximation for quark stars. 
While the GRB data set used is limited nevertheless it 
seems to support the idea that extremely compact objects 
{M oc R^) are behind GRBs activity within our model. 

Acknowledgements. The authors thank J. Schechter, K. 
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